Band- Selective Filter in a Zigzag Graphene Nanoribbon 
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Electric transport of a zigzag graphene nanoribbon through a step-like potential or a barrier 
potential is investigated by using the recursive Green's function method. In the case of the step- like 
potential, we demonstrate numerically that scattering processes obey a selection rule for the band 
indices when the number of zigzag chains is even; the electrons belonging to the "even" ("odd") 
bands are scattered only into the even (odd) bands so that the parity of wavefunctions is preserved. 
The so-called valley- valve effect can be explained by this selection rule. In the case of the barrier 
potential, by tuning the barrier height to be an appropriate value, we show that it can work as the 
"band-selective filter" , which transmits electrons selectively with respect to the indices of the bands 
to which the incident electrons belong. Finally, we suggest that this selection rule can be observed 
in the conductance by applying two barrier potentials. 

PACS numbers: 73.23.-b, 73.63.Nm, 73.40.Gk 



Graphene, a single-layer graphite, is one of the most 
intriguing new materials, which has been studied vigor- 
ously since Novoselov et al. first succeeded in fabrication 
in 2004 [1]. Particularly, its excellent transport proper- 
ties such as the high mobility at room temperature have 
attracted a lot of research interest to apply it in mak- 
ing future electronic devices However, to realize 
graphene-based devices, there is a crucial problem: the 
characteristic gapless band structure makes the control 
of the electric current rather difficult. Many attempts to 
overcome the problem have been done so far, and gave 
rise to various benefical systems such as a graphene quan- 
tum dot [1] , epitaxial graphene on SiC substrates [f| , and 
a graphene nanoribbon [6|. 

We shall focus on a graphene nanoribbon with zigzag 
edges [3] (hereafter we call it "zigzag ribbon" ) . A band 
structure of a zigzag ribbon has two well-separated val- 
leys (conically-shaped curves) K and K' around the ver- 
tices of the first Brillouin zone. In addition, the low- 
energy bands (the lowest conduction band and the high- 
est valence band) are almost flat at the Fermi level due 
to the edge states 0, 0|- Such a peculiar band struc- 
ture has motivated much researchers to investigate the 
electric transport of a zigzag ribbon. Wakabayashi and 
Aoki 0] found that the electric current is almost entirely 
blocked by a barrier potential when the incident energy 
E\ is in the range [0, A] while the barrier height Vq is in 
the range [A, 2A] (2A is the energy separation between 
the top of the next highest valence band and the bottom 
of the next lowest conduction band; see Fig. []Ja)). This 
result is quite different from the case of bulk graphene, 
where incident electrons in the low-energy bands can pass 
a barrier potential of any height (known as the Klein 
Paradox |9(). At first, the polarization of two valleys was 
considered to be the origin of this current blocking effect, 
and then it was named "valley- valve effect" [1(3] by using 
an analogy from the spin- valve effect llj. In the recent 
study [l 21 ] . however, Akhmerov et al. pointed out that 



the origin is not the valley polarization and showed that 
the behavior of the conductance is strongly connected 
with the parity of the number of zigzag chains N. 

This striking current blocking effect has the possibility 
to control the electric current easily. In this letter, we in- 
vestigate the role of a barrier potential applied to a zigzag 
ribbon in more detail, and suggest the way to control the 
current in a new light. The single-orbital tight-binding 
model is employed to describe the electronic states of a 
zigzag ribbon. The Hamiltonian is written as 



(i) 



where the hopping integral ry = r when i and j are 
nearest neighbor sites, and r,j = otherwise. V(yi) is 
the electrostatic potential energy applied at the site i, 
which varies only along the direction of graphene lead 
lines (see Fig. Hfb)). The zero-bias conductance from 
the left lead to the right lead is given by the Landauer 
formula Il3l 
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G = — ^T(/i,z/), T(/i, j/) = |t p 



(2) 



where the summation runs over all incoming (v) and out- 
going (/i) channels, t^ v is a transmission coefficient, and 
T{n, v) represents a transmission probability. We spec- 
ify channels v and /j, by the valley and band indices; the 
right-moving channel in the band n belonging to the val- 
ley K (K') is denoted by K n (K' n ). For example, the 
transmission probability of the scattering process from 
the incoming channel in the band "0" belonging to the 
valley K into the outgoing channel in the band "—2" be- 
longing to the valley K' is denoted by T(K'_ 2 ,K ) (see 
Fig. QJc)). Hereafter we call bands with even (odd) in- 
dices "even" ("odd") bands. For the calculation of the 
transmission coefficients, we adopt the recursive Green's 
function method developed by Ando [3, EH • 
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FIG. 2: (color online) The transmission probabilities for (a) 
v — Ko and (b) v — K\ with N = 30. The insets are the en- 
largements of the low-probability regions. Although we do not 
show T(K^,K ), T(K 4 ,K ), T(K' Al K ) and T{K' 3 , Ki) in the 
diagram, they have a very small yet non-zero value (less than 
~ 10" 4 ), respectively. The others (T(Ki,Ko), T(K 2 ,Ki), 
etc.) are suppressed to 0. 



FIG. 1: (a) The band structure of a zigzag ribbon for the case 
of N = 30. 2 A ~ 3ttt/N is the energy separation between 
the top of the next highest valence band and the bottom of 
the next lowest conduction band, (b) A schematic diagram 
of a zigzag ribbon. The electrostatic potential varies in the 
shaded region. The x axis (y axis) is taken to be perpendicular 
(parallel) to graphene lead lines, (c) A schematic diagram of 
the scattering process from Ko into K'_2 (see text). The bands 
are indexed as shown in the diagram. 



A step-like potential. — Now, we consider the effect of 
a smooth step-like potential described as 

V(y) = V Q(y), (3) 

f (y< -d), 

e(y) = { (1/2) [ S m(ny/2d) + 1) (\y\ < d), (4) 

I i (y > d). 

In the present calculation, we set d = 10a (a is the lattice 
constant) and the incident energy is fixed at E\ — 0.5r 
while Vo is varied from to r. The transmission proba- 
bilities are obtained as a function of e = Ej — Vq. From 
the results of the calculations, we find that scattering pro- 
cesses caused by the step-like potential obey the following 
selection rule: the electrons in the even (odd) bands can 
be scattered only into the even (odd) bands and other 
scattering processes are entirely restricted. For example, 
Figs. [DJa) and [2Kb) show the transmission probabilities 
for v = Kq and v = K\ with = 30. The electrons in 
the band are scattered only into the bands 0, ±2, and 
±4. Similarly, the electrons in the band 1 are scattered 



only into the bands ±1, ±3, and —5. 

The origin of this selection rule is the conservation of 
the parity of wavefunctions [H, 17, 18, [H, |2(J. As is 
shown in Fig. OJa), a zigzag ribbon has the reflection 
symmetry with respect to x = line when N is even. 
Thus, wavefunctions must be either even or odd func- 
tions of x. For example, as is shown in Figs. E^b) and 
[3jc), the wavefunctions of the incoming channel Ko and 
the outgoing channels K^2 and K^4 are even functions, 
while the incoming channel K\ and the outgoing chan- 
nels K'_ 1 and i^_3 are odd functions. Similarly, other 
wavefunctions for even (odd) bands are even (odd) func- 
tions. Since the step-like potential given in Eq. ([3]) can- 
not change the parity of wavefunctions, only the scatter- 
ing processes which preserve the parity of wavefunctions 
can occur. Thus, the transmission probabilities such as 
T(K-2, Ko) have non-zero values and the ones such as 
T(K'_i,K ) are restricted to (see Figs.EJa) andEJb)). 
On the other hand, when N is odd, a zigzag ribbon does 
not have the reflection symmetry and wavefunctions are 
neither even nor odd function of x. Consequently, the 
selection rule does not exist in this case. 

The so-called valley-valve effect [To| can be explained 
by this selection rule as follows. When E\ is in the range 
[0, A] while Vo is in the range [A, 2A], the scattering pro- 
cesses from Ko into K'_ x are required for incident elec- 
trons to penetrate the step-like potential. However, when 
N is even, the wavefunction of Ko is an even function 
while the one of K'_ x is an odd function. Thus, the scat- 
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FIG. 3: (a) A schematic diagram of a zigzag ribbon for even 
N. Open and filled circles represent two inequivalent sub- 
lattices, respectively. One can regard a zigzag ribbon as a 
quasi-one-dimensional system. In this regime, the columns a 
and (3 are treated as quasi sublattices. (b)(c) The wavefunc- 
tions can be written as ^K n (K' ) {x, y) = ipK n (K' ■ We 
show i^K n (K' n ){x) of (b) the incoming channels Ko, Ki on the 
quasi sublattice a in the left lead for E\ — 0.5r, (c) the outgo- 
ing channels K'_ x , K-2, K-z and K-4, on a in the right lead 
for e — — 0.5t. The wavefunctions on the quasi sublattice 
have the same parity as the one of the wavefunctions on a. 
Here N = 30. 



tering processes from Kq into K'_i are entirely restricted 
because it does not preserve the parity of the wavefunc- 
tion, and the current is blocked by the potential. When 
N is odd, Since the selection rule does not exist, the 
current can flow through the potential (l2j |. The similar 



discussion on the valley- valve effect is shown in Ref . [17 . 

barrier potentials. — By calculating the transmission 
probabilities individually, we have just shown that scat- 
tering processes obey the selection rule when N is even. 
However, the existence of the selection rule shown in 
Figs. [2a) and[2jb) cannot be ascertained in experiments 
because the conductance is proportional to the sum of all 
transmission probabilities. Now, we show that the exis- 
tence of the selection rule can be observed in the conduc- 
tance by applying two barrier potentials. 

First, we consider the case of one barrier potential 



V(y) = V {e(y)-&(y-L )]. 



(5) 



To see the role of this potential, we introduce following 
quantities 



P. 



Po = l-Pc 



(6) 
(7) 



where // c (/i ) represents the outgoing channels with even 
(odd) index. The results of the calculations are shown in 
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FIG. 4: (color online) (a) The results of the calculations of 
P c and P . The integral ratios shown in the top (bottom) 
panel represent the values of Ne/N [No/M) (see text), (b) 
Schematic diagrams of band structures and the barrier poten- 
tial described as Eq. ([5]). Here Lo = 50a. The colored regions 
shown in the diagram correspond to the ones shown in (a). 



Fig. BJa). We can see at once that some plateaus ap- 
pear when e > 0. In this region, the contributions of 
the transmission probabilities of the scattering processes 
between the same channels (T(Kq, Kq), T(K%,K\), etc.) 
are major (~ 1) and the others are minor ( < 1CP 3 ) (see 
Fig. Thus, P c and P D are given as P e ~ Af e /Af and 
P a ~ N a /N, respectively. Here N c (Af ) represents the 
number of the right-moving channels with even (odd) in- 
dex in the potential region, and TV = W c + Af Q . 

Let us focus on the behavior in the colored regions 
shown in the diagram. In the red (or light gray) region 
(hereafter called the "even-pass region"), e crosses only 
the band (see Fig. Sib)), and thus P c ~ 1 and P Q ~ 
due to the selection rule [21(. In contrast, in the blue 
(or dark gray) region (hereafter called the "odd-pass re- 
gion"), £ crosses only the band — 1, and thus P Q ~ 1 and 
P c ~ 0. These results mean that a barrier potential in a 
zigzag ribbon plays the role of a "band-selective filter" : 
when e is tuned to be in the even- (odd-) pass region, the 
barrier potential transmits only the electrons in the even 
(odd) bands. 

Next, we append one more barrier potential to the sys- 
tem as follows: 



V(y) 



v [e(y) 

V'[Q{y- 



-O(y-io)] 



Lo-Li-£a)],(8) 



where Vq is fixed so that the first barrier potential be- 
haves as the band-selective filter (i.e. \E\ — Vol < A), 
and V is varied from to t. Introducing two barrier 
potentials like that, one can observe the existence of the 
selection rule in the conductance. For example, when e 
is tuned to be in the even-pass region, the first potential 
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FIG. 5: (color online) (a) Schematic diagrams ol band structures and the barrier potential described as Eq. J5J. (b) The results 
of the calculations of the conductances for Lq — 1/2 = 50a and L\ — 100a. In the upper (lower) panel, Vb is tuned to be 0.426r 
(0.574t) so that e is in the even- (odd-) pass region, and thus the odd- (even-) pass region in the second potential turns to 
the pseudogap region. In the region e' < 0, only interband scatterings occur, which makes the behavior of the conductances 
complicated. 



works as the band-selective filter which transmits only 
the electrons in the even bands. In this case, since the 
incident electrons in the odd bands are reflected by this 
band-selective filter, all the electrons after passing the 
first potential belong to the even bands. Then, one can 
ignore the odd bands in the second potential because the 
scattering processes from even bands into odd bands are 
restricted by the selection rule. This indicates that the 
odd-pass region in the second potential turns to the pseu- 
dogap region (see Fig. EJa); the negligible odd bands are 
shown as dashed lines). In fact, as is shown in the upper 
panel of Fig.[5jb), the conductance is entirely suppressed 
when s' = E\ — V is in the odd-pass (pseudogap) region. 
Similarly, when e is tuned to be in the odd-pass region, 
the even-pass region in the second potential turns to the 
pseudogap region (see the lower panel of Fig.[5^b)). Tun- 
ing the heights of two barrier potentials appropriately as 
above, one can handle the electric current easily, which 
will enables to apply a zigzag ribbon in making electronic 
devices. 

In summary, we have studied electric transport of a 
zigzag ribbon through a step-like potential or a barrier 
potential by using the recursive Green's function method. 
It has been shown that scattering processes in a zigzag 
ribbon obey the selection rule for the band indices when 
the number of zigzag chains N is even. The origin of 
the valley- valve effect is this selection rule. Moreover, 
we have also shown that a barrier potential can play the 
role of the band-selective filter, which transmits only the 
electrons in the bands with either even or odd index de- 
pending on £ = Ei — Vq. Finally, we have suggested that 
the selection rule can be observed in the conductance by 
applying two barrier potentials to a graphene ribbon. 
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